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Recently developed “ray methods" for the asvaptetic uslurion 


of linear partial differential equations are applied to. a system of tira- 


dependent sagneto-ionic equations. These ere equations for E, H, and sha 


° e 
electron current, J. For time-harnonie waves J ean be tliminated and the 


equations reduce te the usual equations of the magneto-Zonie theory. The 


time-dependent equations contain a jarge paremeter A, the plasma frequency, 


which serves as ad expansion parameter. In the “strongly anisotropic” 


case both the gyre frequency and plasma frequency sre large. In the 


singler 


“weakly anisctropic™ case only the plasma frequency is large. 


Both cases can be treated by the general method presented. The weakly 


anisotropic case is analyzed in greater detail. CGne resule of this 


’ analysis is a formula for the rate of rotation of the electric vector 


about. a curved ray in an inhomegeneous medium. The formula reduces to 


that of the well-known Faraday rotation for the special case of a plane 


wave propagating in a homogeneous mediua. 


% 
In recent years aeveptotic getheds haya bean ceveloped fer rhe 


sclution of a larga slaee of problems far Linear partial differenciot 
(equations. ‘thege problems iavolve s pereseter and the mechade prewvide 
One of wore teras of the asyngtetic Sxpanaion, tay fer larze values of : 
the parameter, of the aclution of a given initial or voundary value 
problem. They sre often applicable te problees far wivich ze enac® 
solution method ia known, end ewen for grablems chat cam be solved exs2tly 
it freqiently happens that only the asyastetic expansion of the sclutioa 
‘de sufficientiy ajapie ee be ueeful in practical applications. An 
important thaw ae aeymptotic wethede is characterised by the fact that 
certain curves, or “raya, play a central eehe ia the theory. The rays 
ave {aportant because the functions which make up the various terms of 
the expension satiefy ordinary differcatial squatione elong these curves, 
Often these equations, can be aolved explicitly te yield relatively simple 
formulas chat yieig alisavsete tag inte the physical nature of the 
solutionz. These “ray wethods” are tlosely related to and generalize the 


methods of geometrical optica and Keller's geemmteical theory of diffrac= 


* etom [5]. 


Ta tute pekes we apply ray methods to-the solution of a system 
of partial differential equations that describe sbalvesnaunatie wave 
er ipacetioe dia the ionesphere, By combining Maxwell's equations, in 
their time-dependent form, with a sinple equation of motion for free 
electrons we obtain a hyperbolic system on pakedel differential equations 


fer the electric and magnetic field £ and H, and che eleetron curvest ge 


“> 


wee z a 


in ganeies 2 wa show that fer time-hargende wees we ann eliminires J and 
this syatem reduces to the usual equations of the aaguatedenic *ueor, ih,- 
Therefore we refer to ths equations far 2.8, and J as “time-depeusens 


magnetoienic aguationus". They coutais 4 large parameter, 4, the average 


plassa frequency, which eerves os our expansion parmseter. An equivaicat 


dimensioniese saremater iz i &/e, waers i ig « characteristic length of 
the problem emi c is the speei ef light. tn section 4 we distinguiah tee 
eases, bath of which canbe treated by the general mathod discussed ia 
eacticn 4. In the “ntrongiy snisetrepic” case beth che gyre frequeney end 


plaewa frequency are large peremetera. im the “weekly anisetrosic" cess 


-galy the plagms frequency ts large. 


Tha discussion in section 3 is a beief but self-contained treat- 
ment of the methed developed in (4) and [6] for the auymptotic sulutiea of 
dispersive symaetris hyperbolic ayetens. It includes an improved treatment 
of the colutien of the "transport equations! for the case of 2 double root. 


Thie case te important for the epplication considered here. 


tn section 4 we apply the general resuite of section 3 to the. 

tine-dependent magnetotenie equeticne. Fer simplicity we consider the 
weakly andsotrepde ease and find that thers are three types of nodes. Tha 
“propagating odes" ace the most ieataeaes They contein only “wave 
frequencies” greater than the pleems Scequency, propagate with "group 
velocity" less than c, end are atteauated with “decay exponent" properticnal 
to the ceilision frequency. In addition te the propagating modes there are 
“atanding electric modes” corresponding to an oscillation with the Iecal 


plasms frequency, and stetic magnetic modes. For both of these modes the 


group velocity is sere. ; 
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La Béetion 3 we copadder che preecemtidme: mote amd obi 
aigection of tha field wectera. We iad chet thar ere perpendicular to 
the vey bat estste about it. Gea derive ax ezpidedir forveula for the rate 
of Fotsiien in an Imhesoxyencoug motive, fox which the rave way be curred. 
. Fer the packed CABG CF m time-harmonic place wava propagating ta a 
homogenceus esdiua tha formula reduces to that of che well-keoen Yaraday 


. 


ROCAEAG. 


: Although we have pursued the weakly eniestropia case furcher, ws 


euohadiae that the strongly anisetropic tase is sine included to the general 
_ theory of sectien 3. Vurther anslysie of thie case requires the study of 
the dispersion relation aed cha null vectors ef the “ddoversion mgirix™. 

Tae difiviculties that arise hare are aigebraic anc easentially the same 

as cious thee arige in the usual magueto-tende theacy, where they have tenn 
atalysed in detail {1,2}. It e@ems clear that many of the resulta of the 
steodarcd theery een be applied to ous siparcoeia theery, bus we bave pet 


wndertaken te de thet hare. 


in H.K.S. upita Maxeul Oat] equations for free space take tua fo: 


oxe 2E, sa t s 
® tye) = 4, 
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- Here EB and H are the electric and magnetic field wectore and ¢ and y 
Oe : : * 


the dielectric constant and magnetic permeability of free space. It 
Liows immediately frou (2) that = & (ue) @ ©, hence 2 (3) ie eutomatieslly 
ea fied if it is satisfied et a given ee sey t-0. Thus (3) aay be 
xeted ag a condition on the initial data ot an initial value problem 
T hecunpies field. We shall view (4) es « formula for the 


density p. Then, except for the initial data coadition (3) we need 


r ie consider (1) and (2). 
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In order to deseribe electromagnetic wave propagation in « plasas 


sume that. the ‘current density vector f, consists af two parts, 


feta. |. | (3) 
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motion (feeton’a equatioa) 
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Bere is the slectronic nese amd 9 da (prascribed) external magnetic 


field. We sesume that [Bi<c [R,| asd we have seglected the contribution 


a. af A to the Lorentz force, i.e. to the bracketed cers im (7). ‘The aecand 
* tore tera to {7} te a’ damping force dus te celliaieas, and v ie called 
the collisisn freauenc 1, Thus (6) and (7) yield . 7 


J, = SBE - eu B5 - vl, (8) 
£ B 2 ; 


Wa new intreduce the plasma tre 
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Thus (1), (2) sad (8) becane 
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Tl ese -partial differential equations fer E,H, snd dg wiil be calied the 
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Now (173 yielda 


HF 
in 
mat 


wg = i682 E a 


heace J * ok where so : 


me 


rd “ - on 
oa tcp Hi 122) 
7 
Thus, by eliminating J vom (15-17) we obtain the equation 
eH + ime Ek o sh + Ji, 724} 


YE ~ douk « 6, (263 
which are the usual equations.of the magsete-ienie theory [1,2]. It ia 
interesting to note that whereas (21) isvelvea the relatively simple 
matcix M, in sder to eliminate J wa must introduce the auch mora compiicates 
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[ é “23 #3 
me Se (25) 
Than for any yeoler & * ia exedly seen that 
_, (ZA & xv (26) 
We also intraduse the vector 5 datined by 
go Vest. . 429} 
Then (14) becoses 
dp 7 wUxh + vb > MeoE = 6, (28) 


and tha syetes of equations (123, (13), (28) cae be writcen in the 


bleck matrix uetacies’ 


ae oR) Yeo |izi fa 
a o Hei=lot. (sy 
a 2 a 
[vee 0 gee Ue . 


Each entry in the square matcix on the iefc slide of (29) represents a 
3x3 matrix. Thus, e.g. *e¢ © ved Ty where 13 4a tha 3 « 3 idemedcy metriz. 
Furthermore «ach entry in the colum matrices represents 2 column vecter 
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of dimanaion 3. 
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We sow introduce chrea 9 x 9 matrices a, a, Ao with entries 


: edther sere or one. These matrices can be defined simultanesusly by writing 
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Since 125} igs anti-sypemetric, we ege Crom (350) that fhe watrices A are 


avanctsic and from (31) that A” is positive definite. 3f follows thet our 


differential equations |7). - 


ry 


tm (33) the frequency parametere ¢,y ,Vv and the unit vector J 


~~ 


may be functions of x. It is seavenient ta introduce 2 constant character~ 


rahe mares 


Rene TNS 


istic frequency 4, which may ee defined as a sulteble average of the plasms 
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tn che eeouel we will cascrite a method of ebtaining asymptotic expansion: 
for A+ © of agiutions of hynerbolic systems of the form (74). By 
Satroducing dfimensionices variables, if is not difSisuit ta 6h shad ati 
2 
aguivealent CLerge} cimensionlase parameter ia } wiff’s where ie a 
we 
; ; 3 f i f2, : 5, ote Adicts a 
sharacteriatic length and ce{sp; te the speed of Light. 


Since [eee (343 we have taken both the slaema frequency ¢ and 


the gyre frequency y to be fe 
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take only the plaema frequercy 
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further in sections 4 and 5. 
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We aote that in both cages (35) and (37), 8 ie avnti-symmetric. 


Shig property cf che matrix B 


hich wa shall outiine in the 


is ageential for the succesa of the methed 


next eeetion, {Ita phyeiesi. interpretation 


~ 10 - 


that the large teem AB is: (36) conserves euscgy). Yor Shag re igik 


we Cauhec treet the (dissipative) collision frequency v as a large parmeatar. 
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Asyaptotie theory of hyperbolic systems 


ka {4,6} we considered hyperbelie syetens of the fore 


a 

) 2 

Mu, + Lay, + My + Cy @ fs eo 
Va “vy 


where — Al,...,4° B,C are seooth mug sacrin functions of t end 

x™ (y,.60.8)), plt,x) and £(t,z) are o-dimensional column vectors, ead 

\ is a large pogitive parameter. Here we present 2 byief but seiit-contaiied 
summary of the method of {5,6] for the special case of an asympoticelly 
couservative sysmetric hyperbolic equation, 2.e€- an equation of the form fi} 
for which A? is posicive definite, Al,...,A are Hermitian, and [i i# anti~ 


Hermitian. These cunditiona ere 231 astiefied by the syetem introduced. 


in aection 2. ‘ = 


We consider a formal saymptstic axpansion 


3 
poate (esa) C0) ay (2) A+re@, (2) 


which is te be a solution of (1) with £20. The funetion s is calied the 


phase function and the functions 24] are calied amplitude coefficients. 


We introduce the partial derivatives 


3 
we oe/ge , k= ox, (viel po. 0) (33 
and the dispersion matrix 
. : n = : ; 
Gt ,xw k= J_k A’=4B-0a". (4) 
yeh 


Ry inserting (2) into (i) with f20. and sguutiag ao owes the cuesficlent 


of each power of A, we obtetn the recuraive gyre" -* aqus"2> 


and 


None-trivial solutions are obtaius’ suiy if G te admgeler i.e, if 


the dispersion relation, 


o ahabpaiaciraeEmeTR aH 


det Cedet(J kA’ aiB-wa yo, (7) 
wee] 


is satisfied. Lez 
®eh(t x, i? (8) 
be a root of (7) of multiplicity q. This means that for weh there exist 


4 ilnearly independent null vectors r such that 


A 
f 
Gr,* Q, cd Phys «34 e £9} 
me ; = 
’ 5 oem 4 3 a 
Since, for real k, the matrix })  & 4 ~i8 te Hermitian, and A” is positive 
vel 


definite, there exist real roots 1), the aum of their multiplicities q 
le equal to m, er d the combined eet cf & Bull vacters ig Linaarly tedepencent. 
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The vectors yican be so chasen that 


on 


o 3 x 
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Equation (8), considered ae a first order partial differential 
equation for a(t,x) [eee (3}i can be solved by rhe mathod of characteristics 


{3}. i£€ we introduce the group velocity vector & with components 


g teh / ak . 7 7 wel, ee ay mn : (11} 


& , 
If a and Bare column vectors the inner product (3,5) is defined by 
pt 
(a,b) = ) a4 by. The etary denotes the complen conjugate. For Matrices 
a hese ; 


it denotes the Hermitian conjugate, 
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(18) 


From (63 we see that co is in the null gpace ef G, hence has 


an expansion of the form 
P| 
goes softy . (39) 
qm 


n (5) we may ingers (9) for go Tf we then take the inner 


belo 
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produc: of the resulting equatioa with ro, the left sits vanishes end the 
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Sometimes the term “ray” is alec used for the projections of theae curves 
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Thus we tice that the coelficienta 91,..4,0, in the expension (19) ‘cor 
i sCe 
the Les noeiitcient g ef (21; satiefy the system of q Firat 
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r L1eey <3 eae are ar as saney 6% % 
er 4 ciiferentisi equations (transport eguations) 124% 


Sho higher order amplitude casfficients 2 


se 


, (433) are determined by simila~ 


quetiogs whisk ace derived in 14}. 


In order te discuss methods of solving the system ox; transport 
equations (24) we must firet introduce the “ray transformation.” In any 


given syoblen we will have an n~parameter family ¢f soluticus 


« xft.y), &™ k(t, 7) of the characteristic equatione (12). aad hence 
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Here we have fatreduced the expansion ef the determinant by cc faccors of 


no ate ‘ 
the 2 row. Since a determinant vanishes when two Powe are Ene sane, wa 


have the identity, 


ye cot 2x4 o 36 (25) 
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By differentiating (25) and using (26) avd (iz) we obtain 


eee a A a*n4 CES ee a ee 22,5 8 en dx4{ j 
3 dt” Lay te cof oy, § aD Ge? | ar, oF ay, | 


wo 8x co “2 
aw To (Xtye § cs + 27 
"sex, Ge? Tru x og + (27) 


This identity which relates the jacobian of the ray transformation and 


the divergence of the group velocity vector will be used shortly. 


We sow impose che following condition om the coefficient 
matrices of (1): 
(C + ox = a ~ Eg ot i yea? = x6 ga 3 B,GL,...4q. (28) 

“vee 

The condition states that the qxq matrix whose entries are given by the left 
side of (98) reduces to a acalar. The canteen is trivially satiefied if 
g=i.' For the aeopeten sass equations of section 2, the matrices AY are 
constant hence (28) reduces te aconditier on the matrix €, which we shall 


examine in section 4. Now from (22), (173, (27) and (28) we obtain 
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Thus we have obtained a formula for the energy density at any point t on 
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@ vay pe xlt,y) | in terme of ite valve eat a fixed seint t, on the same rsy. 
The cepencence of ail quanticies in (32} om the parameters y which “label” 


the rave ig not explicitly indicated. 


The energy density is an important scalar function that characterizes 


the leading amplitude coefficient z, of the asyiptetic solution (2). in seme 
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esses we can obtain furcher informaticu abot 2 by integrating the traaspor: 


equations (24). The moat important cases are rhosa for which the multiplicity 


q ig one or two. 
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Thea P is anti-hermitian. 
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